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1 INTRODUCTION 



ABSTRACT 

We have acquired deep, H-band, imaging for a sample of 286 Virgo cluster galaxies 
with Br ^ 16 mag and extracted surface photometry from optical g,r,i,z Sloan Dig- 
ital Sky Survey images of 742 Virgo Cluster Catalog galaxies, including those with 
H-band images. We confirm the detection of a dip in the luminosity function indica- 
tive of a discontinuity in the cluster galaxy population; the dip is more pronounced at 
redder wavelengths. We find, in agreement with earlier works of TuUy & Verheijen and 
ours for Ursa Major cluster galaxies, a clear dichotomy between high and low surface 
brightness galaxy disks. The difference between the low and high brightness peaks 
of Virgo disk galaxies is '^2 H-mag arcsec"^, significantly larger than any systematic 
errors. The high surface brightness disk galaxies have two distinct classes of high and 
low concentration bulges, while low surface brightness galaxies have only low con- 
centration bulges. Early- type galaxies exhibit a similar structural bimodality though 
offset from that of the spiral galaxies towards higher surface brightnesses. Both the 
early- and late-type structural bimodalities are uncorrelated with colour or any other 
structural parameter except, possibly, circular velocity. Random realizations of realis- 
tic surface brightness profiles suggest that a bimodal distribution of effective surface 
brightness is unexpected based on normal distributions of bulge and disk parameters. 
Rather, the structural bimodality may be linked to dynamical properties of galaxies. 
Low angular momentum systems may collapse to form dynamically important disks 
with high surface brightness, while high angular momentum systems would end up as 
low surface brightness galaxies dominated by the dark halo at all radii. The confir- 
mation of structural bimodality for gas-rich and gas-poor galaxies in the high-density 
Virgo cluster as well as the low-density UMa cluster suggests that this phenomenon 
is independent of environment. 



The range of surface density profiles for galaxies in a cluster 
is a telltale of its dynamical history. For instance, an excess 
of cuspy profiles over, say, the mean field galaxy distribution, 
is indicative of recent merger activity. Conversely, extended 
exponential profiles are representative of quiescent evolu- 
tion over long timescales (e.g., Toth & Ostriker 1992). The 
study of galaxy light profiles has a rich history with seminal 
early contributions from de Vaucouleurs (1948; 1959), Sersic 
(1968), and Freeman (1970). De Vaucouleurs (1948) estab- 
lished that the surface brightness (hereafter SB) profiles of 
early type stellar systems may follow a strongly concen- 
trated distribution now referred to as the "de Vaucouleurs" 
or "r^''*" profile. De Vaucouleurs (1959) later reported the 
ubiquity of exponential SB profiles in disk galaxies while 
Sersic (1968) proposed a generalized fitting function that 
encompasses the exponential and de Vaucouleurs profiles. 



Based on photographic images for 36 disk and SO galax- 
ies. Freeman (1970) postulated that the distribution of disk 
central surface brightness (hereafter CSB) , defined as the in- 
tercept of an exponential disk fit at r = 0, peaks at /io=21.65 
B mag arcsec"^. Deep, wide-field CCD, galaxy surveys have 
since revealed a rather continuous distributions of CSBs 
from high surface brightness (HSB) to low surface brightness 
(LSB) galaxies (de Jong & Lacey 200ofl- 

Simulations and 

observations of galaxy structure have thus far suggested a 
continuous range of properties over, rather than fundamen- 
tal departures from, the HSB and LSB regimes. However, 
based on near-infrared (NIR) observations of the Ursa Ma- 
jor (UMa) cluster, TuUy & Verheijen (1997; hereafter TV97) 



^ Besides the small number statistics. Freeman's study was biased 
by the low threshold of photographic plates to bright galaxies. 
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inferred that the distribution of disk central surface bright- 
nesses could be bimodal. 

Based on a sample of 62 Ursa Major cluster galaxies 
observed at optical and NIR bands, TV97 argued that the 
distribution of disk central surface brightnesses, ^o, was 
not continuous as expected from simple structure forma- 
tion models, but rather bimodal. The bimodality of the 
/io distribution was not convincingly observed at optical 
(BRI) bands, which TV97 argued might be due to ex- 
tinction effects at optical bands. At Kq TV97 could as- 
sert the existence of two peaks in the SB distribution at 
17.28 K' mag arcsec"^ for HSB galaxies (corresponding to 
"Freeman's law") and at 19.69 K' mag arcsec"^ in the LSB 
regime. The surface brightness of the gap was centered at 
~18.5 K' mag arcsec"'^. Inclination and extinction correc- 
tions to the optical surface brightness profiles would bring 
a CSB bimodality to light in these bands as well. TV97 
stressed, however, that the SB bimodality at optical bands 
could plausibly be an artifact of their extinction correction 
which applies only to HSB galaxies. Still, the fact that the 
K'-band is nearly insensitive to dust and that the distribu- 
tion of /xo,K' values is clearly bimodal is a strong case for 
TV97's argument of structural bimodality in spiral galaxies. 

TV97 divided their sample into galaxies with and with- 
out a significant near-neighbour, to identify an environmen- 
tal effect. The small UMa membership however means that 
the isolated and crowded sub-samples contain only roughly 
30 galaxies each. Despite the statistical limitations, TV97 re- 
ported evidence for an enhancement of the bimodality in the 
isolated sub-sample. The majority of galaxies with interme- 
diate /xo measurements were found to have near neighbours. 
This led TV97 to suggest the existence of two stable dy- 
namical configurations of baryonic and dark matter, likely 
induced by the environment, leading to the distinct LSB and 
HSB populations. 

Comparing light distributions with total mass distri- 
butions inferred from HI synthesis maps allowed TV97 to 
further classify these two states. In LSB galaxies, dark mat- 
ter dominates the potential at all radii; gas and stars have 
sufficient angular momentum that they have not settled to 
the core and are orbiting in response to the dark matter dis- 
tribution. In HSB galaxies, baryonic matter has dissipated 
its energy and transferred angular momentum sufficiently 
to become self-gravitating over the central ~ 2 exponential 
disk scale lengths (Courteau & Rix 1999; Button et al. 2007, 
hereafter D07). A gap between LSB and HSB galaxies may 
suggest that galaxies avoid a situation where baryonic and 
dark matter have comparable dynamical influence in the in- 
ner disk. The higher incidence of intermediate surface bright- 
ness (hereafter ISB) galaxies with near neighbors suggests 
that tidal influences might drive a galaxy from the LSB to 
the HSB state. 

Although TV97's study provided several valid argu- 
ments for a structural bimodality of cluster galaxies, it still 
suffered a few shortcomings. First, the UMa cluster mem- 
bership is small. With only 62 galaxies, a structural study 



2 The K' filter was described by Wainscoat & Cowie (1992). It 
resembles the 2MASS Ks, with the primary difference being that 
it excludes the longest wavelength part of the K atmospheric win- 
dow in order to reduce thermal emission. 



is plagued by statistical uncertainty (Bell & de Blok 2000). 
Second, TV97 estimated the disk central surface brightness, 
/io, by fitting an exponential profile to the outer disk in 
order to avoid the bulge component. The uncertainty in 
this "marking-the-disk" technique depends on the subjec- 
tive interpretation of the bulge size and the disk fit baseline. 
This caveat was however examined carefully by McDonald 
et al. (2008) and found not to be the cause of the observed 
bimodality. Finally, TV97's result could be unique to the 
UMa cluster and not representative of the true nature of 
disk galaxies in clusters or the field (de Jong & Lacey 2000). 
Because the UMa cluster lacks early type galaxies, the ob- 
served bimodality could simply be due to a missing morpho- 
logical class of galaxies. It is, however, our impression that 
none of these objections are stronger or more convincing 
than the claim of bimodality itself, as we shall verify below 
for Virgo cluster galaxies. 

Following the release of the Sloan Digital Sky Sur- 
vey (York et al. 2000; hereafter SDSS), there have been 
many reports of bimodality in the distributions of galaxy 
color, star formation and clustering properties (Strateva 
et al. 2001; Blanton et al. 2003; Kauffmann et al. 2003; 
Baldry et al. 2004; Brinchmann et al. 2004; Balogh 
et al. 2004) and even reports of trimodality in galaxy concen- 
trations (Bailin & Harris 2008). Colors, star formation and 
clustering of galaxies are all intimately linked with a tran- 
sition baryonic (light-weighted) mass of 3xlO^"M0 (Kauff- 
mann et al. 2003). Red galaxies, at the upper end of the mass 
scale typically have low current star formation rates and 
are found primarily in cluster or high-density environments, 
while the less massive blue galaxies have high star formation 
rates and are found primarily in the field (Dressier 1980). 
These results have catalyzed several new variants of galaxy 
formation models and evolution. For instance, an appealing 
explanation for any star formation bimodality involves two 
phases of gas accretion in galaxies via hot and cold modes 
(Birnboim & Dekel 2003; Dekel & Birnboim 2006). The hot 
mode is the standard picture of gas that is shock heated to 
the virial temperature during gravitational collapse of the 
gas and halo. This shock-heated gas eventually cools and 
settles into a galaxy (Rees & Ostriker 1977; White & Rees 
1978; White & Frenk 1991). The cold mode consists of cold 
filaments penetrating far inside the halo (Birnboim & Dekel 
2003). These can occur in low- mass halos where virial shocks 
are absent and cold gas is accreted quasi-spherically. For ha- 
los with larger masses, the gas is shock heated to the virial 
temperature everywhere except where cold, dense gas fil- 
aments penetrate the virial radius. As the mass increases, 
only filaments of increasingly dense cold gas can penetrate 
the halo without being shock-heated. Keres et al. (2005) 
used semi-analytic models to predict that the transition halo 
mass where hot accretion becomes more efficient than cold 
accretion is ^--^3 x 10^^ M©. This concords with the transi- 
tion baryonic mass, 3 x 10^° M0, in the color/SFR bimodal- 
ity, assuming a universal baryonic to dark matter ratio of 
~10% (Zavala et al. 2003). Keres et al. (2005) also found 
that low-mass halos in dense environments have enhanced 
hot accretion. The cold mode would thus apply mostly to 
late-type, disk galaxies while the hot mode pertains mostly 
to early-type, spheroidal, galaxies. This simplistic picture 
appears to match observations well, but realistic star for- 
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mation prescriptions are needed for more robust data-model 
comparisons. 

1.1 Plan of action 

We now wish to test TV97's claim of structural bimodal- 
ity in a different environment and using a larger sample. 
Thanks to its proximity, large size, and broad morpholog- 
ical coverage, the Virgo cluster is the next logical cosmic 
structure to consider as we revisit TV97's analysis. If con- 
firmed, we can also assess whether the reported bimodality 
of surface brightnesses, for spiral galaxies at least, is in any 
way related to the SDSS color bimodality discussed above 
and whether the global trimodality of galaxy concentrations 
for all Hubble types (Bailin & Harris 2008) is also observed. 

The organisation of this paper is as follows: In !j2l 
we discuss briefly our database for the 742 Virgo Clus- 
ter Catalog (Binggeli et al. 1985; hereafter VCC) that are 
found in the SDSS 6th Data Release (Adelman-McCarthy 
et al. 2008; hereafter DR6). Out of these, we have defined 
a complete, magnitude-limited, sample of 286 VCC galaxies 
with Bt ^ 16 for which deep H-band photometry was ob- 
tained. The full description of our data set, and details about 
the data acquisition and reduction methods are given in Mc- 
Donald et al. (2009; hereafter "data paper"). We introduce 
in 131 specific parametric and non-parametric quantities that 
will permit a comprehensive understanding of the structure 
of Virgo cluster galaxies of all types. In Appendix A, we 
determine the effects of various errors on the distributions 
of various non-parametric quantities in ^X] We present final 
results in §31 and examine their possible interpretations in 

m 

We assume in this paper a distance to all Virgo cluster 
galaxies of 16.5 Mpc or m-M=31.18 (Mei et al. 2007). At 
that distance, 1" — 80 pc. 



2 VIRGO SAMPLE 

The study of the unbiased distribution of galaxy surface 
brightness requires volume completeness. A solution to this 
issue is to study galaxies at a common distance in a clus- 
ter. Our sample is drawn from the Virgo Cluster catalogue 
(Binggeh et al. 1985; hereafter VCC). The VCC catalogue 
contains 2096 galaxies within an area of ^^140 deg^ on the 
sky, centered on the galaxy M87 at a~12''25™ and 5~13°. 
The VCC is asserted to be complete down to a limiting ab- 
solute magnitude of Mb ~ — 13 and to contain many objects 
as faint as AIb ~ — 11 (Binggeli et al. 1985). 

Our first task was to acquire deep, optical, photom- 
etry for as many VCC galaxies as could be found in the 
SDSS/DR6. A further spatial cut was made following Tren- 
tham & TuUy (2002) to reject 40% of the cluster that is con- 
taminated by the W, W' and M background groups. The W, 
W' and M groups were identified by de Vaucouleurs (1961) 
and Ftaclas et al. (1984), respectively. W and M lie about 
twice Eis far away as Virgo and W' lies about 50% further 
than Virgo. These and related structures in the Virgo South- 
ern Extension all lie in a flattened plane close to the su- 
pergalactic equator, contaminating the western edge of the 
Virgo Cluster. Trentham & TuUy (2002) discuss this prob- 
lem and describe how the projected area of the cluster can 
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Figure 1. Two-dimensional map of the Virgo cluster of galaxies. 
Colored points are galaxies that belong to our H-band sample 
and that lie at a distance of ~ 16.5 Mpc. Grey points are galaxies 
in our VCC/SDSS sample, while black points are the remaining 
galaxies in the VCC. The size of the points scales with the total 
luminosity. The concentric dashed circles are projected galacto- 
centric distances (1,2 and 3 Mpc) around the center of the Virgo 
cluster at M87. The dashed red circle corresponds to a distance 
of 6° from M87. 



be separated between zones with negligible and serious con- 
tamination. 

Images from the SDSS/DR6 were extracted for a total 
of 742 VCC galaxies obeying these spatial cuts. We will refer 
to this sub-sample as the "VCC/SDSS" sample. 

In order to ensure that our Virgo catalog reaches be- 
low the LSB peak of UMa (TV97), we must achieve com- 
pleteness at least down to Mb = —16.45 mag (assuming a 
distance of 15.5 Mpc to UMa). In order to ensure that no 
intermediate surface brightness galaxies are missed, we can 
define a complete sub-sample that includes all VCC galax- 
ies with AIb ^ —15.15 mag (i.e. Bt ^ 16). We also make 
an additional two spatial cuts to ensure that we only in- 
clude bound cluster members: 1) all galaxies identified as 
"background" by Binggeli et al. (1985) are removed and 2) 
all galaxies further than 6° from the center of the clus- 
ter, defined as the position of M87, are removed. These 
magnitude and spatial cuts leave us with a total of 303 
VCC galaxies centered on M87. A further 8 VCC galax- 
ies (1068,1217,1258,1355,1665,1768,1889 and 2096) with re- 
cessional velocity measurements, Vrad > 3000 kms~^, and 
the galaxies VCC0723 and VCC0991, with significant fore- 
ground stars, were also excluded. The following seven VCC 
galaxies (530,950,1052,1287,1571,1822 and 1992) could also 
not be detected at H-band (see ^2.2|l , and were left out of the 
sample (this has no effect on our final conclusions). We are 
left with a final, complete, sub-sample of 286 VCC galaxies 
with Bt ^ 16 , covering a wide range of luminosities and 
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Figure 2. Distribution of galaxy morphologies in the full VCC catalog (black), in the VCC/SDSS sample (grey), and in the H-band 
sub-sample (red). The percent ratios are the number of objects in a given morphological bin for each of the two samples divided by that 
in the full VCC. 



morphologies. For reasons that wiU soon be clear, we will 
refer to this complete sub-sample as our "H-band" sample. 

Fig. [1] shows the distribution of all VCC galaxies in 
black, as well as the 742 VCC/SDSS galaxies in grey, and 
the H-band sample in multi-colors. Fig. [5] shows the distri- 
bution of galaxy morphologies, as taken from the NASA Ex- 
tragalactic Database (hereafter NED). The broad morpho- 
logical coverage is important to ensure that no distribution 
of galaxy structural parameter is biased by morphological 
segregation. Fig. [2] shows that Virgo is intrinsically rich in 
dE, SO and Im galaxies (see also Mei et al. 2007). This is 
not due to a bias in our sample but simply to the nature of 
the Virgo cluster (as represented by the VCC). 

The gas-rich VCC galaxies are represented at approxi- 
mately the same levels in our VCC/SDSS and H-band sam- 
ples relative to the full VCC with completeness between 11% 
and 42%. The bright magnitude limit of the H-band sample 
also implies a higher number of early-type galaxies relative 
to the later types. However, we will show in !j4] that the 
bimodality of surface brightnesses is detected in each mor- 
phological bin (early or late type) and that it is thus unlikely 
due to the morphological make-up of this specific cluster or 
our sampling of it. 

Following TV97's study of the UMa cluster, we have 
sought to obtain optical and especially near-IR photome- 



try as well as dynamical measurements for all of our H- 
band sample galaxies. The SDSS ugriz imaging for this sam- 
ple will yield a distribution of luminosities, surface bright- 
nesses, scale lengths and concentrations, as a function of 
optical wavelength. The optical colors will enable a compar- 
ison of any surface brightness bimodality, if present, with 
the observed SDSS galaxy color bimodality (e.g., Strateva 
et al. 2001) and the ability to determine if one is simply 
a consequence of the other. The near-IR images are how- 
ever essential to uncover the true distribution of galaxy sur- 
face brightnesses, luminosities, and other galaxy parameters 
largely free of extinction by dust. We describe the extraction 
of the optical and NIR imaging data below. 

While most critical, our collection of dynamical parame- 
ters for VCC galaxies is still in progress and will be reported 
elsewhere. 

2.1 SDSS Photometry 

We have extracted calibrated ugriz images from the 
SDSS/DR6 for 742 VCC galaxies, including the 286 galax- 
ies in our "H-band" sample. Surface brightness profiles and 
total luminosities were obtained for these galaxies in all five 
SDSS bands by, first, performing isophotal ellipse fitting to 
the i-band images according to the methods of Courteau 
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(1996), and then, applying the i-band isophotal solutions 
to the images in the other SDSS bands. The latter ensures 
that color gradients extracted from all SDSS images are 
computed from the same matching isophotes (MacArthur 
et al. 2003). u-band images were consistently shallower than 
the griz bands and were thus discarded. Sky levels for back- 
ground subtraction and the photometric zero-points for cah- 
bration were obtained from the SDSS image headers and the 
SDSS archives, respectively. The remainder of the profile ex- 
traction technique is identical to that used for the near-IR 
photometry, as described below. 



Table 1. Summary of H-band observations for the 286 Virgo 
cluster galaxies in our sample. 



Tel - Camera 


Targets 


Collected 


Avg Seeing 


UH88" - ULBCAM 


52 


04/2005 


1.2 ± 0.2 


UH88" - ULBCAM 


16 


04/2006 


1.5 ± 0.3 


UH88" - ULBCAM 


31 


04/2007 


1.3 ± 0.2 


UH88" - ULBCAM 


23 


03/2008 


1.0 ± 0.2 


UKIRT - WFCAM 


31 


07/2008 


1.1 ± 0.1 


CFHT - WIRCAM 


34 


02-06/2008 


1.1 ± 0.1 


GOLDMine 


79 




2.2 ± 0.9 


2MASS 


20 




2.6 ± 0.1 



2.2 NIR Image Collection 

Due to practical constraints, new, deep, H-band imaging 
that would sample well below the putative surface bright- 
ness bimodality scale of TV97 could only be obtained for 
a smaller sample of VCC galaxies. This is the magnitude- 
limited "H-band" sample of 286 VCC galaxies with SDSS 
imaging described above. 

Deep H-band imaging for some VCC galaxies is already 
available from the Two Micron All-Sky Survey (Skrutskie 
et al. 2006; hereafter 2MAS^ and from the GOLDMineQ 
database (Gavazzi et al. 2003). We were able to secure H- 
band imaging for the remainder of the H-band sample with 
the detectors ULBCAM at the UH 2.2-m telescope, WF- 
CAM at UKIRT and WIRCAM at CFHT over the period 
2005-2008. Some deep K-band imaging was also available for 
a few galaxies in 2MASS and GOLDMine. 

H-band images from GOLDMine were kindly provided 
by G. Gavazzi. Calibrated 2MASS galaxy images were ex- 
tracted from the online database. Many of the 2MASS and 
GOLDMine images were not deep enough for our purposes. 
Whilst adequate for large HSB galaxies, the high 2MASS 
brightness threshold (typically fin = 21 mag arcsec"^. Bell 
et al. 2003; Courteau et al. 2007; Kirby et al. 2008) hmits the 
use of those data bases for deep extragalactic studies. Like- 
wise, just a handful of GOLDMine images were deep enough 
to properly separate the bulge and disk light. We have de- 
fined the relative depth criterion, Q, as the ratio of the maxi- 
mum extent of the H-band surface brightness profile divided 
by that of the SDSS i-band profile: Q = rmax.Nm/rmax.i, 
where rmax is the radius where the surface brightness error 
exceeds 0.15 mag arcsec"'^ (see McDonald et al. 2008 for de- 
tails). Wherever possible, we impose Q > 0.75 for our NIR 
data. Ultimately, 20 2MASS and 79 GOLDMine galaxy pro- 
files were deemed useable for our study. A final 187 Virgo 
cluster galaxies needed new observations. These new and 
existing observations are summarized in Table [1] 

The surface brightness profiles for the deep GOLDMine 
and 2MASS images were measured using the same tech- 
niques as our new NIR images to ensure uniformity for the 
entire database. Further details regarding the data reduc- 
tion process and quality are presented in McDonald et al. 
(2008; 2009). 



http : //www . Ipac . caltech.edu/2niass/releases/allsky/ 
http://goldjnine.mib.infn.it/ 



3 SURFACE BRIGHTNESS PROFILE 
ANALYSIS 

Surface brightness profiles were extracted for the 742 
VCC/SDSS galaxies at griz bands and for the 286 H-band 
galaxies. In this section, we examine the various parametric 
and non-parametric properties derives from these profiles. 

McDonald et al. (2008) already considered the paramet- 
ric (i.e. model dependent) analysis of UMa spiral galaxies in 
order to test, and ultimately confirm, TV97's claim of bi- 
modality. The current, larger, Virgo cluster sample however 
includes galaxies of all morphologies. While a parametric 
approach to model the shape of complex light profiles in- 
volves a multi-component decomposition (e.g., MacArthur 
et al. 2003; McDonald et al. 2008), a non-parametric ap- 
proach is free of model assumptions and reveals difi'erent as- 
pects of galaxy structure. We explore both approaches below 
and it will be shown later that surface brightness bimodal- 
ity is essentially independent of the method of light profile 
analysis. 



3.1 Parametric Quantities 

Spheroidal and flattened galaxy systems have traditionally 
been modeled as the sum of a bulge and disk components 
(see MacArthur et al. 2003; McDonald et al. 2008 for more 
details). The ID light profile of a galaxy disk is typically 
parametrized as an exponential function: 



Idir) = Joexp {--^l , 



(1) 



where 7o and h are the disk central surface brightness 
and scale length, respectively. Meanwhile, the projected ID 
bulge light profile is better modeled as a Sersic function 
(Sersic 1968): 



Ib{r) = 7e exp < ~bn 



l/r. 



1 



(2) 



where Te is the half-light radius, le is the surface bright- 
ness at that radius, and n is the Sersic shape parameter. 
With n = 1, the Sersic function reduces to the exponential 
function. 

In addition to the bulge and disk, we consider other 
components that may affect the light profile such as com- 
pact nuclei, spiral arms and disk truncations. Ignoring these 
may result in large errors in the bulge and disk parameters 
(McDonald et al. 2008). 
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B/D decompositions were performed only for the H- 
band and optical griz light profiles of the 286 galaxies in 
the "If-band" sample (McDonald et al. 2009). 



3.2 Non-Parametric Quantities 

Non-parametric quantities are measured directly from the 
surface brightness profile, with no prejudice for any assumed 
model. The only assumptions inherent to the non-parametric 
measurements below are: (i) that the total galaxy light is an 
extrapolation of the light profile to infinity, and (ii) that 
inclination estimates and extinction corrections are valid. 
Non-parametric quantities allow a direct, unbiased compar- 
ison of galaxies across the full Hubble sequence. 



3.2.1 Concentration, C28 

The galaxy light concentration is a measure of the rela- 
tive light fraction between the inner and outer parts of 
the galaxy. Unlike the B/D ratio, which relies on a model 
for the light distribution, C28 is a straightforward, model- 
independent morphological indicator. The concentration, 
C28 is defined as (e.g., Kent 1985; Courteau 1996): 



C28 



V r2o 



(3) 



where r8o and r2o are the radii enclosing 80 and 20 percent 
of the total light. For a pure exponential function (n = 1 in 
Eq. ©), C28=2.8. 

Concentration indices are a function of wavelength and 
while fractional radii depend on inclination and extinction 
corrections, we expect the ratio r8o/»'20 to be roughly in- 
dependent of projection effects. We justify this assumption 
below. 



3.2.2 Fractional Radii and Surface Brightness 

Fractional parameters refer to specific quantities measured 
at radii that contain specific fractions of the total galaxy 
light. The effective radius, re, is the radius that encloses 
half of the total light. While re was already introduced in 
the context of a parametric profile function (Eq. [2}, it is 
formally defined in a non-parametric way as: 



I[r)2TTrdr 



■I 



I{r)2TTrdr. 



(4) 



The effective surface brightness, Je, is the surface brightness 
at re; le = ^('■e) (Me = — 2.51og/e). Unlike the extrapolated 
central surface brightness of the disk, Iq, the effective sur- 
face brightness le is non-parametric, making it ideal for the 
comparison of mean surface brightness levels for galaxies of 
varying morphology; le (or fie) applies to all galaxy types 
rather than 7o (or ^0) which is restricted to spiral disks 
(Jo = 1.678/e for pure disks). Our discussion about the vari- 
ations in galaxy surface brightness profiles will indeed rely 
on re and He. 

Another fiducial light marker is the fractional radius 
defined as: 




Figure 3. Distribution of H-band effective radius, re,H, and 
concentration, C28 for the Virgo H-band sample as a function 
of i-band axial ratio (top), and against similar i-band mea- 
surements in the middle section. The point /line types repre- 
sent: dE-dS (red, open/dashed), E-SO (red, closed/solid), Sa-Sb 
(green, closed/solid), Sc-Sd (blue, closed/solid) and Irr (blue, 
open/dashed). The black, dotted, line in the middle windows is 
the slope unity correlation. The bottom windows are population 
histograms. 



where Tj, contains some percentage (x) of the total light. 
fix is the surface brightness in mag arcsec"^ at r^;. A final, 
non-parametric, measure of surface brightness is the average 
surface brightness interior to some radius rx, (Ix), defined 
as: 



JJ^^ I{r)2TTrdr 



(6) 



/ I(r)2TTrdr = 

/o 100 



I{r)2TTrdr 



(5) 



{fi)x is the magnitude equivalent of {I}x. 

Fig.|3]shows the dependence of the H-band re^H and C28 
on the j-band axial ratio, b/a for the H-band sample. There 
appears to be a correlation between re,H and b/a, but it is 
driven largely by the spheroidal nature (high axial ratios) 
of compact galaxies (small re). Any correlation between re 
and b/a is significantly weakened if compact galaxies are 
excluded. The insensitivity of C28 to projection effects (top 
right window) is even clearer. The respective corrections for 
projection on r2o and rgo roughly cancel out. We also see in 
Fig. O that there is a good linear correlation between i- and 
H-band scale radii and concentrations. 

Ultimately, concentrations, fractional radii and surface 
brightnesses were computed for all 742 VCC/SDSS at gri 
wavelengths and for all 286 "H-band" galaxies at H-band. 



4 RESULTS 

We now examine the distribution of structural parameters 
of Virgo cluster galaxies, starting with an analysis of the 
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gas-rich systems in the H-band sample for comparison with 
TV97's similar analysis. We then derive the distribution of 
effective surface brightnesses and the optical and NIR lumi- 
nosity functions for the full H-band and VCC/SDSS sam- 
ples, with all morphologies considered. 

We will find for the Virgo gas-rich galaxies in the H- 
band sample, a bimodal distribution of the extrapolated disk 
central surface brightness, /io,H, and thus bolster the simi- 
lar claim for UMa galaxies by TV97. For the complete sam- 
ple of gas-rich and gas-poor Virgo galaxies, we determine 
the distribution of efi'ective surface brightness, fie (indepen- 
dent of B/D decompositions) and find compelling evidence 
for brightness bimodality in each morphological group. We 
present in i|4. li the distributions of fio and, in i]4.2l various 
fractional and average surface brightness measures. We com- 
pute in i|4.3l the optical and NIR luminosity functions and 
compare these to the recently published optical luminosity 
function of Virgo (Rines and Geller, 2008) and field NIR 
luminosity function from UKIDSS (Smith et al. 2008). In 
i|4.4l we present the distributions of various scale radii and 
the galaxy concentration, C28. The latter is compared to the 
distribution of concentrations for SDSS galaxies by Bailin 
and Harris (2008). In i|4.5l we examine the bivariate distri- 
butions for most of the structural parameters addressed in 
this section. 



4.1 Disk Central Surface Brightness 

In order to compare with the study of UMa cluster galax- 
ies by TV97 and McDonald et al. 2008, we restrict our 
H-band sample to disk galaxies. B/D decompositions, as 
in McDonald et al. (2008), were performed for this sub- 
sample of 166 VCC disk galaxies. The B/D fits, which in- 
clude optional bulge, nucleus and spiral arm components, 
were applied to each griz and H-band SB profiles. The de- 
rived /xo values were corrected for projection effects to their 
face-on value in the absence of extinction. Thus, we write: 
/io = Mo — 2.5Ca log(fe/a), where b/a is the measured ax- 
ial ratio of the outermost isophote at a given band in each 
galaxy. We take Ca = 1 for the dust transparent case. This is 
likely a fair assumption at H-band, the focus of our analysis, 
but less adequate for the griz bands. We keep our analysis 
free of dust correction for now. 

The distribution of inferred fig values is shown in Fig. U 
There is a clear dearth of galaxies at /Iq u ~19 mag arcsec"^ . 
This result for later-type VCC galaxies matches very well 
that determined by TV97 and ourselves for UMa disk galax- 
ies with a minimum in the number of galaxies at /Xq ~ 
18.5 mag arcsec^^. Similarly, we find excellent agreement 
with TV97 in the location of the HSB and LSB peaks at 
Pq^^=17.85±0.15 mag arcsec"^ and /ioj^=20.27±0.4 mag 
arcsec"^, respectively, with a peak separation of A/ig ^^=2.4 
mag arcsec"^. Considering a typical H-K~0.2 for galaxy 
disks, the distributions for Hq^h in the Virgo and UMa clus- 
ters are very similar. 

Using a statistical F-test, we can reject the hypothe- 
sis of a normal distribution for /Xg f{ in favor of a bimodal 
distribution with a confidence level of 85%. Fig. [5] shows 
which galaxies contribute mostly to different surface bright- 
ness levels. As expected, early-type disk galaxies dominate 
the HSB peak, while late-type disk galaxies and irregulars 
are present in both the HSB and LSB peaks. 




22 20 M 

/Xq^ [mag arcsec^^] 

Figure 4. Distribution of inclination-corrected disk central sur- 
face brightnesses, fig ^ for 166 VCC spiral and irregular galaxies 
(hatched). The red line shows the results from the analyses of 
UMa galaxies by TV97 and ourselves (McDonald et al. 2008), 
while the green line shows the two-Gaussian fit to the VCC data. 
The F-test value demonstrates that a single-Gaussian fit can be 
rejected with a confidence of 85%. The transition between the two 
brightness peaks is at /Xq ~ 19mag arcsec"^. 
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Figure 5. Inclination-corrected disk central surface brightnesses, 
fi'g with morphological type for 166 VCC galaxies. Filled and 
open black squares represent HSB and LSB galaxies, respectively. 
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Figure 6. Distribution of inclination-corrected disk central sur- 
face brightnesses, fi^, from optical (SDSS) to near-IR bands for 
166 spiral and irregular Virgo cluster galaxies from the H-band 
sample. The red histogram in the i-band panel refers to the 272 
disk galaxies in the VCC/SDSS sample. The numbers on the right 
are the statistical F-test probability for a 2-component Gaussian. 

Let us now examine the distribution of inclination- 
corrected /io with wavelength. Fig. |6] shows the different 
distributions from g (top) to K (bottom). The gap between 
the two SB peaks grows as a function of wavelength since the 
HSB peak brightens at a faster rate with wavelength than 
the LSB peak. This is explained by the higher dust content 
in HSB galaxies, in contrast with their relatively transparent 
(C^ = 1) LSB counterparts (TV97), as well as the relative 
colors of the constituent stars (HSB galaxies tend to be red- 
der than LSB galaxies). Dust obscuration is less effective at 
longer wavelengths and, as a result, we see more deeply along 
any given line-of-sight thus increasing the observed surface 
brightness. Recall that we have not applied any correction 
for obscuration (Ca=1) to surface brightnesses and the bi- 
modality disappears as we consider shorter wavelengths. It 
is possible that the correct choice of Ca would restore the 
bimodality at optical wavelengths, as suggested by TV97, 
but any such complications can be avoided by restricting 
our analysis to NIR wavelengths only, as we do here. 

We can ask if the bimodality in the disk CSB is corre- 
lated with any bulge structure. Fig. [7] shows the same dis- 
tribution as in Fig. |4] but now with a second /io,£r distribu- 
tion including only galaxies with small bulges. We define a 
"small" bulge as one which contributes < 25% of the total 
galaxy light. Indeed, we see that if we remove those galaxies 
with bulges which contribute more than 25% of the total 
luminosity, bimodality is preserved. Though the F-test con- 
fidence level drops, this is primarily due to the reduction in 
sample size - the ISB gap is very well preserved. 

Fig. [S] shows the relationship between fio and disk scale 
length, h, at H-band. At intermediate luminosity, there can 
be both HSB and LSB galaxies, depending on the disk sur- 



Figure 7. Histogram of /ip ^ for 166 spiral and irregular VCC 
galaxies. The red, hatched histogram includes only those galaxies 
with very little contribution from the bulge component (bulge-to- 
total ratio, B/T, < 25%). This shows that bimodality is not the 
result of bulge-dominated vs bulgeless systems. 



20 I I I I 




Figure 8. Distribution of fio as a function of disk scale length, h, 
for the sub-sample of 166 spiral and irregular VCC galaxies. Red 
dashed lines represent constant luminosity; our detection limit is 
~1O^L0 at H-band. Filled and empty squares correspond to HSB 
and LSB galaxies respectively, as in Fig. [5] 
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Figure 9. Correlation between the parametric brightness, fio, 
and non-parametric brightnesses fie and (fJ.)^ for the 166 gas-rich 
VCC galaxies. The point types are as follows: open red circle - 
SO, green circles - Sa & Sb, blue circles - Sc & Sd, open blue 
circles - irregulars. The dashed lines are fig = — 1.822 and 
Mo = ~ 1-124, representing a pure exponential disk. 

face density (faint/high no & large h or bright/low fig & 
short h). The fact that two galaxies with the same lumi- 
nosity can have wildly different surface brightnesses leads to 
the belief that some mechanism likely related to the initial 
halo angular momentum (Dalcanton et al. 1997) prevents 
LSB systems from collapsing to the same densities as the 
HSBs. This figure makes clear that for a given luminosity the 
LSB galaxy must be more extended than the HSB galaxy 
and that no LSB galaxy exceeds 10^^ L©. The distribution 
of disk scale lengths also shows some evidence of multiple 
peaks, though the statistical significance is low (68%), none 
of which are correlated with the structure in ^o.^f. 

We have shown thus far that the NIR central surface 
brightnesses, /Xq ^ of disk galaxies in the Virgo and UMa 
clusters are distributed bimodally. Our B/D decompositions 
confirm a result that is already well-known - that galaxy 
bulges come in two types: cuspy (high concentration) and 
cored (low concentration). HSB disks harbour both types 
of bulges, while LSB disks only harbour low-concentration 
bulges. 

We now expand our analysis to take full advantage of 
the diverse Virgo cluster population by considering early- 
type and dwarf galaxies as well. 

4.2 Effective Surface Brightness 

Unlike the UMa cluster of galaxies, Virgo is rich in giant and 
dwarf early-type galaxies. The mean density of Virgo is ~5 
times that of UMa. Many of the VCC galaxies have no ob- 
servable disk component and therefore our analysis of galaxy 
properties based on /lo would be moot. We consider instead 
two more versatile quantities: fie, the galaxy effective surface 



Figure 10. fj.o — fJ-e and fio — {fJ.)e residuals as a function of 
concentration. 

brightness and (fJ.)^, the mean surface brightness, defined 
as the average surface brightness within re (§4.2.2). These 
non-parametric quantities can be measured for all galax- 
ies. Fig. 1^ shows the comparison for the 166 gas-rich VCC 
galaxies of with fie and (fi)^. There is a strong correla- 
tion between fio and fie with the expected zero-point offset 
{fie = /io + 1.822) for pure exponential disks). The scat- 
ter increases for earlier-type galaxies. This scatter is further 
examined in Fig llOl where we confirm that the differences 
between fio,fie, and (fi.)^ are a function of morphology or, for 
simplicity, concentration. That is, for higher concentration 
the scatter in the fio — fie and fio — {fi}e relations increases. 

We have also considered the average surface brightness 
within re, {y)ei a-s a measure of a galaxy's characteristic 
surface brightness. Use of this parameter is motivated by a 
study of late- type field galaxies by de Jong & Lacey (2000; 
hereafter DLOO). DLOO studied the distribution of {fi)e for 
1000 Sb-Sdm field and cluster galaxies which, they reported, 
is not bimodal. The right side of Fig. [5] shows that {fi) ^ cor- 
relates well with fiQ , though with slightly larger scatter than 
fie- From this point on, we will adopt fie as our standard 
measure of surface brightness for three reasons: 1) For spiral 
galaxies, fie scales directly with fio, 2) fie can be measured 
for any galaxy morphology, and 3) fie is independent of any 
assumptions about the shapes of the galaxy bulges (if any) 
and disks. We also adopt a geometric correction for fie in 
all galaxies: fil = fie — 2.5\og{b/a)r^, where {b/a)r^ is the 
axial ratio at the effective radius. Regardless of what surface 
brightness measure we use, and how we correct for projec- 
tion effects, we will show in §6 that the overall shape of the 
distribution of surface brightnesses is preserved. 

We show in Fig. [11] the distribution of fil with wave- 
length for all 286 H-band VCC galaxies. The differential ef- 
fects of extinction with wavelength are well known (though 
poorly understood) and likely the cause for the smoother 
contributions at shorter wavelengths. The fie distribution of 
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Figure 11. Distribution of cfTcctivc surface brightness, fi],, in 
5 bands for all 286 VCC galaxies in our sample. The red his- 
togram in the i-band window is the distribution of for all 742 
VCC/SDSS galaxies; that distribution is strongly bimodal with 
an F-test of 95%. Troughs are significant in the brightness distri- 
butions from i to H bands. 



the nearly dust-free LSB galaxies should not change drasti- 
cally from g to H, modulo a color term due to stellar popu- 
lations. The HSB peak from (j to H will, however, be altered 
by effective dust obscuration. The faint tail of the HSB peak 
will thus be stretched to fainter values, effectively filling any 
intrinsic trough between the HSB and LSB peaks (as shown 
by TV97 for UMa galaxies). For this reason, we shall now 
rely solely on NIR surface brightness measurements for the 
remainder of our analysis. 

Fig. [12] shows the distribution of /i^ ^ measurements 
for our sample of 286 H-band VCC galaxies. Characteristic 
structure can be seen in most of the morphological bins, as 
well as in the sample as a whole. The separate distributions 
of E and SO galaxies are each weakly bimodal. The com- 
bination of these morphological types (top right panel of 
Fig. I12|) gives a stronger bimodality for h].^h- The bimodal- 
ity of for gas-rich galaxies emerges clearly from the 

HSB and LSB peaks for the early-type (Sa-Sbc) and irreg- 
ular galaxies; the distribution of Sc-Sd galaxies shows the 
least features and is least abundant of all VCC types. This 
results in a strong bimodality for all the gas rich galaxies 
(middle-right panel of Fig. \12\ . If we compare the distribu- 
tions of fJ.\^H for gas-poor and gas-rich galaxies, the lower 
surface brightness spheroids line up with the HSB peak for 
late- type and dwarf galaxies. This results in an apparent tri- 
modality of /le for the complete H-band sample (lower right 
panel of Fig. I12p . However, we stress that the trimodality 
is in fact a superposition of an early- and late-type SB bi- 
modality. The distribution of (/i)|,£f in Fig. 1131 shows similar 
trends. However, Fig. [TS] illustrates the sensitive nature of 
binning for this sort of argument. For instance, the distri- 
bution of {fi)l,H for Sc-Sd is here roughly bimodal when it 
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Figure 12. Distribution of inclination corrected effective surface 
brightness, /le , for 286 VCC galaxies of varying morphology. The 
red vertical dotted lines correspond to the locations of the ESB 
(extreme surface brightness), HSB and LSB peaks. The upper- 
right panel shows all gas-poor galaxies, the middle-right panel 
is the sum of all the gas-rich galaxies, and the bottom-right is 
the sum of all morphological. The sum of the two bimodal distri- 
butions for the gas-poor and gas-rich galaxy types respectively, 
results in a trimodal distribution of all the VCC galaxies in our 
sample. The numbers in the upper and middle right panels refer 
to the F-tost confidence for bimodality. 
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Figure 13. Distribution of average surface brightness, {fi)^, for 
galaxies of varying morphology. The numbers in the upper and 
middle right panels refer to the F-test confidence for bimodality. 
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was flat in Fig. 1121 This illustrates the need for large sam- 
ples. We are still reassured that the distributions for or 
(/i)e for the sum of both the gas-rich and gas-poor galaxies 
show very similar (bimodal) features. Figs. ll2)ll3l show that 
a structural bimodality is present in both early and late- 
type galaxies with a tail towards the brightest galaxies. A 
statistical F-test confirms the rejection of a single peak in 
favor of two distinct peaks with 71% and 82% confidence 
for early- and late-type galaxies, respectively. We can con- 
sider another significance test. In order to fill the early and 
late- type brightness gaps, an additional 15 and 29 interme- 
diate SB galaxies would be needed, respectively. However, 
the current number of galaxies in these troughs is 24 and 
39, respectively, so the required additions to achieve fiat SB 
distributions would correspond to 3 and 4.5 standard de- 
viations. This argument, based only on Poisson statistics, 
complements nicely the F-test (which assumes a Gaussian 
shape) . 

FiEs. ll2)ll3l tell an exciting story; not only is a structural 
bimodality present in spiral galaxies, but also in spheroids. 
Equally intriguing is that spheroids of lower surface bright- 
ness have identical surface brightnesses as HSB spirals, pos- 
sibly pointing to a common formation and/or structural 
mechanisms. In order to better understand this result, we 
now examine the correlation of surface brightness with vari- 
ous other parameters and what role it plays in fundamental 
structural relations. But first, let us introduce another un- 
expected result, this time in the NIR luminosity function, 
to shed more light on the structural bimodality in both gas- 
poor and gas-rich galaxies. 
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Figure 14. Optical and NIR luminosity functions for our sample 
of 286 Virgo cluster galaxies. In the r-band panel, the open trian- 
gles correspond to Petrosian magnitudes taken directly from the 
SDSS for comparison with Rines & Geller (2008). In the i-band 
panel, the closed triangles correspond to magnitudes computed 
for all of the 742 VCC/SDSS galaxies. Open circles refer to mag- 
nitude bins suffering from incompleteness at the low-luminosity 
end. 



4.3 Near-IR and Optical Luminosity Functions 

McDonald et al. (2008) showed that the K-band luminosity 
function for Ursa Major has a noticeable dip at intermediate 
luminosities. We now examine this claim for Virgo galaxies 
with both optical and near-IR luminosities. Fig. ll4l shows the 
luminosity functions of our 286 VCC galaxies in the g,r,i,z 
and H bands. There appears to be a dip in the luminosity 
function at optical wavelengths, which develops into a broad 
trough at H-band. This trough clearly separates two distinct 
luminosity classes. 

Rines and Geller (2008) obtained the same result for 
their study of the faint end of the Virgo LF (see black tri- 
angles in r-band window of Fig. I14p . but dismissed it on 
grounds that SDSS data may be fiawed at luminosities where 
the dip is observed. It would however be challenging to think 
of such a systematic bias in the SDSS pipeline reductions. 
Furthermore, our luminosities are measured from SDSS im- 
ages but fully independent of any other SDSS data product 
(e.g., Petrosian luminosities and radijfl). The dip in the Virgo 
luminosity function is also observed at H-band which further 
confirms that it is not an artifact of the SDSS pipeline. Our 
luminosity function and that of Rines & Geller are a striking 
confirmation of each other, establishing the reality of the dip 
in the luminosity function of Virgo cluster galaxies. 

The dip in the NIR luminosity is not reproduced in the 
UKIDSS luminosity function for field galaxies (Smith et al. 

® Rines and Geller (2008) used "model" luminosities from the 
SDSS pipeline. 



2008). This could indeed mean that there is no trough in 
the luminosity function of field galaxies and that the lack of 
intermediate-luminosity galaxies in UMa and Virgo is purely 
a cluster-driven phenomenon. However, the UKIDSS lumi- 
nosity function suffers incompleteness at luminosities fainter 
than Mk = —20, while the low-luminosity peak in the UMa 
luminosity function is at Mk = —19.5. This fact, along 
with the difficulty of achieving volume completeness in field 
galaxy samples, could conspire to hide any putative trough 
in the NIR luminosity function if it is, in fact, independent 
of environment. 

We will see in §4.5 and Fig.[221that the high- luminosity 
peak is comprised of "extreme" surface brightness (hereafter 
ESB) galaxies and the low-luminosity peak is comprised of 
LSB galaxies. However, the HSB galaxies are distributed 
bimodally in luminosity, implying a bimodality in sizes for 
these galaxies since they all have nearly the same surface 
brightness. The combination of Figs 1121 and 1141 shows that 
the majority of the low-luminosity HSB galaxies are compact 
dE galaxies, while the majority of the high- luminosity HSB 
galaxies correspond to traditional HSB spirals. 

Finally, Fig.[T2]shows the multi-dimensional variation of 
various structural parameters with optical and NIR luminos- 
ity for the full sample of 286 VCC galaxies. The dumbbell- 
shaped distribution of points in all plots of luminosity versus 
He or (/i)e shows clearly the separate galaxy populations. 
The L-fi relation appears to be very similar whether we use 
He or (/i)e, although the scatter is somewhat reduced at fie- 
Finally, there is a strong evolution of most of these scaling 
relations with wavelength that will be studied elsewhere. 
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Figure 15. Various structural parameters in X plotted against total absolute luminosity in Y for four optical bands (griz) and the H-band 
for 286 VCC galaxies. Structural parameters from left to right are as follows: He d - disk effective surface brightness, - Scrsic index 
of the bulge, log - half-light radius, fig - effective surface brightness, - average surface brightness. Histograms of the luminosity 

function are shown in the far right panels. The red lines represent the binned median using the same binning as the histogram. 



4.4 Concentration 

In Appendix A, we have attempted to predict the observed 
distribution of C28 based on normal distributions of fid, rd, 
fib, rb and n. In Fig. [16] we show the distribution of mea- 
sured concentrations for the 286 VCC galaxies binned by 
morphology as well as for the complete H-band sample. For 
comparison, we show the predicted distribution from our 
models (Appendix A). 

The predicted and observed distributions are remarlc- 
ably similar, though the observation of a second peak at 
high concentration with C28 > 4.0 was not predicted by 
our simulations of isolated galaxies. The overall distribution 
peaks at C28 ~ 3, corresponding to pure exponential disks. 
The distribution terminates rather sharply at C28 ~ 2.0, 



suggesting a lower limit for how centrally concentrated the 
baryons in a galaxy can be. In order to understand the high- 
C28 peak, let us look at the distribution of C28 for different 
morphological types. Fig. 1161 shows that the Sc-Sd, Irr and 
dE galaxies populate a single peak at C28 < 4, while E, SO 
and Sa-Sbc galaxies occupy two distinct peaks near C28 ~ 3 
and C28 = 4.5. Thus, there is a concentration 

bimodality only in early-type (bulge-dominated) galax- 
ies. The difference between the model and observed distri- 
butions is due in large part to the fact that our model dis- 
tribution of bulge surface brightness is normal, while the 
one we observe is bimodal. The transition at C28 = 4 repre- 
sents the transition from LSB to HSB bulges. We combine 
Figs.[T2]and[T6]in Fig.[l7]to examine whether the early-type 
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Figure 16. Distribution of concentration, C2S, for 286 VCC 
galaxies of different morphologies. The red histogram in the Sa- 
Sdjlrr window refers to the distribution of C28 for UMa galaxies 
from McDonald et al. (2008) for galaxies of the same morphol- 
ogy. The red histogram in the bottom-right window is based on 
100,000 model surface brightness profiles, discussed in §5.4. 
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Figure 17. Correlation between H-band concentration, C281 ^^nd 
effective surface brightness, /i* ^. The dashed vertical line is an 
arbitrary transition at C28 = 3.8 from low- to high-concentration, 
while the two dashed horizontal lines delineate the LSB, HSB and 
ESB classes. 



bimodality in concentration corresponds to the bimodality 
in surface brightness. 

We see again in Fig. [T7] that the two HSB peaks are 
comprised primarily of E, SO and Sa galaxies, while the LSB 
peak is comprised primarily of dwarf, irregular and late- 
type spiral galaxies. Most importantly, we see two isolated 
clouds: one containing high-concentration HSB and ESB 
galaxies and the other containing low-concentration HSB 
and LSB galaxies. Again, LSB galaxies never have high con- 
centration and ESB galaxies never have low concentration. 
This provides some insight into the structural bimodality 
in spheroids. It is possible that galaxies with C28 > 4 can 
only form through major mergers, which naturally explains 
the lack of both LSB and star-forming galaxies in the high 
concentration peak. We can also verify that the trimodality 
observed in SDSS concentrations by Bailin & Harris (2008) 
is not detected here. 

4.5 Bivariate Distributions 

In order to assess if the brightness bimodality extends to, 
or even stems from, other galaxy parameters, we now con- 
struct a number of bivariate distributions for [s,^ vs other 
fundamental structural parameters such as color, luminos- 
ity and dynamical mass. 

The color-surface brightness bivariate distribution is 
certainly of great interest in light of the color bimodality 
of SDSS galaxies (e.g., Strateva et al. 2001). It is of obvious 
interest to ask if the transition in galaxy colors matches the 
transition in galaxy surface brightnesses. A careful examina- 
tion of Fig [18] defuses this idea. The fact that the blue and 
red peaks each include all three SB peaks demonstrates that 



the structural bimodality is indeed different from the color 
bimodality. The same distribution but for g-H colors (Fig. 
I19p reveals a different color bimodality (resulting from stel- 
lar population effects that will be investigated elsewhere). 
The fact that the red peak alone is comprised of all three 
SB classes again confirms that surface brightness trends and 
color bimodality are independent phenomena. We summa- 
rize our findings in Fig. [20] with the distributions of gr — i, 
(7— H and \i\ n for individual morphologies. There is a strong 
correlation between color and surface brightness for elliptical 
galaxies (and weaker for SOs), with redder ellipticals having 
higher surface brightness than blue ellipticals. There appears 
to be no correlation between color and surface brightness for 
the other types, suggesting again that the surface brightness 
bimodality is independent of stellar populations. 

The correlation of /i^ with rotational velocity gives 
greater hope of establishing a link with the observed bi- 
modality in /Xe. Indeed, we note in Fig.[5T]hints of bimodality 
for Virgo galaxies in the distribution of the maximum galaxy 
rotational velocity, Vmai. With high-quality, heterogeneous 
(combination of HI linewidths and Ha rotation curves) dy- 
namical information available from HyperLEDA (Paturel et 
al. 2003) for only 102 of the 286 VCC galaxies, the observed 
weak correlation (r=0.53) between Vmax and /ie.w may not 
be surprising. A homogeneous collection of resolved rotation 
curves for the complete sample of galaxies will be needed to 
establish a dynamical extension of the galaxy structural bi- 
modality. In a similar vein, TV97 found for a subsample 
of 35 UMa galaxies that the surface brightness correlated 
strongly with the slope of the rotation curve, a subject we 
address in Appendix 2. Yet again, a much larger sample is 
needed for a conclusive analysis. 

Finally, we observe the distribution of luminosities vs 
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Figure 18. Bivariate distribution of H-band surface brightness 
and g—i color for 286 VCC galaxies. Point colors and types are the 
same as for Fig. 1171 Dotted lines represent individual Gaussian 
fits while dashed lines are the sum of these fits. 
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Figure 19. Bivariate distribution of H-band surface brightness 
and g-H color for 286 VCC galaxies. Point colors and types are the 
same as for Fig. 1171 Dotted lines represent individual Gaussian 
fits while dashed lines are the sum of these fits. 



Figure 20. Distribution of surface brightness and colors for 286 
VCC galaxies for different morphologies. 
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Figure 21. Bivariate distribution of H-band surface brightness 
and fully corrected Ymax for 102 VCC spiral galaxies. The Gaus- 
sian fits in the upper panel correspond to the full distribution of 
286 VCC galaxies. Dotted lines represent individual Gaussian fits 
while dashed lines are the sum of these fits. 



u in Fig. 1221 There is a clear trend between luminos- 
ity and surface brightness, as expected. Although there is 
scatter in both variables, any bimodality in log Vmax , as ob- 
served in Fig. [21] ought to be reflected in log L by virtue 
of the velocity-luminosity relation of galaxies. Furthermore, 
since ~ (A')e' ''■'^'^ (/^)e ~ -^'i/2/7r'"e we would expect the 



bimodality to be reflected in the distribution of luminosi- 
ties, unless the SB bimodality is due to the relation Ljr^ . 
Fig. 1221 shows a mild bimodality in L, but which is uncorre- 
lated with the stronger bimodality in /ie. 

We have analyzed bivariate distributions involving 
color, rotation velocity, luminosity and surface brightness. 
From Figs. I18II22I we flnd no strong correlation between the 
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Figure 22. Bivariate distribution of H-band surface brightness 
and total H-band luminosity for the 286 Virgo cluster galaxies 
with H-band photometry using a distance to the Virgo cluster of 
16.5 Mpc. Dotted lines represent individual Gaussian fits while 
dashed lines are the sum of these fits. 

structure seen in the distribution of surface brightnesses and 
that of any other structural parameter. There is no correla- 
tion between surface brightness and color, with both HSB 
and LSB galaxies occupying the red and blue peaks. Figs [211 
& IBll suggest that the source of the bimodality in surface 
brightness may be linked to galaxy dynamics and not stellar 
populations. TV97 also proposed a dynamical explanation 
for the structural bimodality. In the following section we of- 
fer a tentative explanation for this exciting result in a more 
global context. 



5 DISCUSSION 

In §4, we have verified that the distribution of surface bright- 
nesses is bimodal in UMa and Virgo cluster galaxies. Before 
we discuss possible explanations for this effect, let us in- 
vestigate any possible biases in the data measurements and 
calibrations that could conspire to produce the observed bi- 
modalities. Environment, which we discuss below, may also 
play a role in the structural bimodality. Various angles to 
address any cause for bimodality - field vs cluster galaxies 
or optical vs IR data, or both - will then be explored. We 
conclude with plausible explanations for the observed galaxy 
surface brightness bimodality. 

5.1 An Artificial Bimodality? 

We here investigate any potential source of bias from our 
measurement and calibration techniques. In our data pa- 
per (McDonald et al. 2009), we quote estimates for the 
error induced by our photometric calibration of ~0.1 H- 
mag arcsec"^. In Figs. [4^ and 1121 we observed that the 




Figure 23. Distribution of /i* ^ for 286 NIR VCC galaxies 
with an over-subtraction of the sky of 0.05% (red), and under- 
subtraction of 0.05% (blue) and the nominal sky subtraction 
(black). 

separation between the LSB and HSB peaks is ~2 H-mag 
arcsec"^. In order to artificially cause this gap, a surface 
brightness error greater than 0.5 H-mag arcsec"'^ on both 
LSB and HSB peaks must be invoked. This error would also 
be systematic in order to "evacuate" ISB galaxies; bright 
LSB galaxies would have to be made fainter and/or faint 
HSB galaxies would have to be made brighter. We will ex- 
amine below three possible processes that could induce such 
a large error: inclination correction, sky subtraction, and 
surface brightness measurements. We will attempt to show 
that, although far from negligible, these errors can be con- 
trolled and should not bias the shape of surface brightness 
distributions. Finally, we will examine the possibility that 
either of the SB bimodalities are due to regions with clus- 
tered HSB or LSB galaxies (i.e. an infalling group). 

5.1.1 Sky Subtraction Errors 

In Appendix A, we demonstrated how a sky measurement 
error of ±0.01% can result in an artificial truncation or up- 
turn in a deep surface brightness profile. However, while LSB 
galaxies are highly sensitive to sky errors, a systematic over- 
or under- subtraction of the sky could hardly enhance the 
ISB gap. In Fig. 1231 we show the distribution of fie if the 
images for all the sample galaxies are systematically over- 
or under-subtracted by 0.05% (a huge error at H-band). It 
is clear that, the net effect is to shift all galaxies towards 
LSB, not to create a gap. 

5.1.2 Inclination Corrections 

In §4.1, we discussed deprojection corrections to measured 
surface brightness quantities. At NIR wavelengths, galaxies 
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Table 2. Location and separations of the three SB peaks for SB 
measurements. 
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Figure 24. Distribution of effective surface brightness for 128 
late-type VCC galaxies using different inclination corrections. 



are nearly transparent and the correction for deprojection is 
purely geometric. 

Since is measured from the disk fit, the geometric 
correction is trivial, provided the inclination of the disk is 
known. However, the quantity /le, which depends on both 
bulge and disk hght, is less easily corrected. We thus make 
use of two different methods to correct for projection. 
First, we can write the corrected effective surface brightness 
of the whole galaxy, fj,\, as: 



— 2.5 log 



(7) 



where (b/a)^ is the axial ratio at the effective radius, r^. 
This correction (Method A) is straightforward, since we 
know the ellipticity as a function of radius from isophotal 
fits (§3.2.3). 

The second inclination correction (Method B) requires 
the decomposition of the SB profile into bulge and disk com- 
ponents and deprojection of the two components separately. 
Knowing the inclination of the disk, we can deproject the 
disk fit (assuming full transparency). We further assume 
that the bulge requires no deprojection. This yields the cor- 
rection for the global nl: 



-2.5 log 



Ibulgeire) + Idiakire) ( 



(8) 



where (b/a)^ is the axial ratio of the outermost isophote. 
The effective radius is measured from the total hght profile 
and the contributions from the individual bulge and disk 
profiles at this radius are used. This correction is intuitively 
appealing as the bulge should require less of a deprojection 
correction than the disk. The assumption of a spherical bulge 
may however be inadequate, especially in the presence of 
a bar. However, we find that the isophotes in the central 
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regions of most of our galaxies that exhibit bulges are, in 
fact, fairly round. 

Fig. [24] shows the result of applying the two afore- 
mentioned inclination corrections. We find very little dif- 
ference between methods A & B. As expected, the cor- 
rected values to face-on orientation for /i^.i/ Stie systemati- 
cally fainter. Most importantly, regardless of the correction 
used, bimodality is always preserved. As there is little differ- 
ence in the resulting distributions, we prefer Method A for 
its simplicity and empirical, rather than parametric, quality. 
Furthermore, Method B is less desirable since it may yield 
large errors in cases when the disk fit is unreliable or if the 
assumption of spherical bulge is unwarranted. 



5.1.3 Surface Brightness Measurements 

Ideally, we wish to avoid any sensitivity to the bulge hght 
as well as to localized radial structure. Such a parameter for 
pure disks is fj,o, but we do not have an analog for both disks 
and spheroids. Indeed, He is affected by local structure while 
(h)^ encompasses bulge light. Of interest to us is to investi- 
gate if the use of a local vs global estimator of surface bright- 
ness at different galactocentric radii should cause bimodal- 
ity. We can compute fix and {fJ,)x ^ in §3.2.2. At short radii, 
we expect both hx and (/i)x to be biased by any bulge or 
bar, while at large radii fix should suffer mostly from SB 
fluctuations and sky subtraction errors. Fig. 1251 and Table [2] 
show the results of our exercise for all 286 VCC galaxies in 
our sample at H-band. The result is interesting: regardless 
of the specific measure of surface brightness, whether local 
or global, computed at small or large radii, the three SB 
peaks are always preserved. 

As can be seen from Table [21 the separations between 
the ESB and HSB peaks, and the HSB and LSB peaks are 
roughly equal. The remarkable result that the peak sepa- 
rations (ESB-HSB) and (HSB-LSB) are nearly identical at 
most SB levels suggests that there may be a single mecha- 
nism responsible for all three peaks. 

Thus, we have shown that both the early- and late- 
type structural bimodality persist regardless of the surface 
brightness measurement, inclination correction, and sky sub- 
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Figure 25. Surface brightness distributions for local and global measurements at different galactocentric radii, fix refers to the surface 
brightness at the radius containing X% of the total light, while (fJ,)x refers to the average surface brightness interior to that radius. The 
blue dashed lines are multiple Gaussian fits that correspond to the superimposed bimodal distributions of early- and late-type galaxies. 
The red solid line is the sum of the Gaussians. 



traction error. Let us then address possible causes for the 
observed galaxy bimodalities, convinced that they are not 
an artifact of our methodology. 



5.1.4 Spatial Surface Brightness Distributions 

It is conceivable that one of the three surface brightness 
peaks could be due to the infall of a small group of galax- 
ies with a narrow surface brightness distribution. In order 
to investigate this issue we turn to Fig. 1261 which shows 
the spatial distribution of galaxies in the four different SB 
peaks: HSB & LSB early-type galaxies, and HSB & LSB 
late-type galaxies. It is clear from Fig.[26]that both the HSB 
and LSB spiral and irregular galaxies are spread uniformly 



throughout the cluster. Likewise, most of the spheroids are 
distributed uniformly throughout the cluster, with an over- 
abundance of HSB spheroids near the cluster center. Thus, 
although there is a constant infall of galaxies into the Virgo 
cluster, it does not appear that this is causing an artificial 
peak in the surface brightness distribution. 



© 2009 RAS, MNRAS 000.[llf27l 



18 M. McDonald, S. Courteau & R. B. Tully 



20 



15 



u 

0) 
Q 



10 



5 



15 



u 

0) 

Q 



10 



5 



I I I I I I I I I I I I I 



... •>••.•.•..- 

. ' • • 

* *. •* Late — type 



I I I I I I I I I I 



• • • • 

•• .. . J,-.-. . • •;!: 



•I 



• *. •/. Early — type 
* HSB 



J I I L 



I I I I I I I I I I I I I I I I I I 
• % 

M. * ■ 

• •• V--. • 



• *• •..*••? Late — type 

^' HSB 

I I I I I I I I I I I I I I I 




•I 



• *. •* Early-type 



• ESB 

I I L_ 



J L 



190 185 
R.A. (°) 



190 185 
R.A. (°) 



Figure 26. Spatial distribution of VCC galaxies in four surface brightness bins: LSB late-type (/^^ ^ > 21 mag arcsec HSB late-type 
(l^e H ^ '^^ ™3'g arcsec"^), HSB early-type (^i* ^ > 18.5 mag arcsec"^), and ESB early-type ^ < 18.5 mag arcsec"^). The smaller 
gray points represent the full sample of 286 Virgo cluster galaxies. 



5.2 Cluster vs Field Distributions 

Freeman's law (1970) in its modern form is understood as 
the maximum surface brightness any galaxy disk can 
but it does not preclude the existence of fainter systems. The 
general picture today is that the surface brightness function 
of galaxies has a tail at low surface brightnesses that extends 
to arbitrarily low values (de Jong & Lacey 2000). Previous 
studies that have attempted to quantify the distribution of 
galaxy surface brightnesses have focused mostly on field sys- 
tems (e.g., McGaugh & de Blok 1995; de Jong & Lacey 2000; 



° Freeman's original sample was heavily biased towards high sur- 
face brightness systems, hence the narrow reported range of sur- 
face brightnesses. 



Driver et al. 2005) . The effect of environment on galaxy sur- 
face brightnesses, if any, would thus be missed. 

In their pioneering study, TV97 observed that the struc- 
tural bimodality is strengthened if their UMa sample is re- 
stricted to isolated galaxies. For UMa, an isolated galaxy was 
defined as a galaxy that has not encountered another galaxy 
during the lifetime of the cluster. Since UMa is dynamically 
young, this definition is a reasonable one. However, since the 
crossing time of the Virgo cluster is less than one-tenth of a 
Hubble time (Trentham & Tully 2002), it is unreasonable to 
perform a similar analysis for Virgo cluster galaxies. Every 
galaxy in the cluster should have had at least one significant 
encounter within the last Gyr, which means that we cannot 
divide the sample into isolated and non-isolated galaxies. 

Various studies have reported that LSB galaxies form 
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primarily in isolated environments (e.g., Rosenbaum & Bo- 
mans 2004) likely from l-2a peaks in the initial density fluc- 
tuation spectrum. However, these authors found that local 
LSBs are almost exclusively found in clusters and filaments. 
If LSB galaxies form from l-2a peaks, they would migrate 
from voids to high-density regions. This does not explain 
how dynamically young galaxy clusters, such as UMa, could 
have a high abundance of LSB galaxies. Instead, perhaps the 
cluster environment enhances LSB formation which would 
explain the current location of the LSB population as well 
as the SB bimodality seen in UMa and Virgo. 

For the UMa and Virgo clusters, there is a bimodality 
of disk surface brightness in spiral galaxies as well as bi- 
modality in effective surface brightness for early-type galax- 
ies. Since previous studies of field galaxies have not shown 
any evidence for bimodal surface brightness distributions, we 
might therefore think of structural bimodality as a cluster 
effect. However, it should be noted that most of the studies 
of the surface brightness distribution of field galaxies have 
focused strictly on optical imaging for lack of available com- 
plete similar data bases at IR wavelengths. We showed in 
Figs. [6] & [TT] that the bimodality in disk and total surface 
brightnesses is enhanced in the NIR, due to uncertain and 
significant corrections for dust at optical wavelengths. None 
of the major NIR data bases available today (e.g., de Jong 
1996; Jarrett et al. 2003; Grauer et al. 2003; MacArthur 
et al. 2003; Smith et al. 2008) cover the full range of surface 
brightness to sample the LSB peak. For this reason, it is 
unclear whether the structural bimodality observed at NIR 
wavelengths is indeed a cluster effect since it has only been 
observed in galaxy clusters thus far, or whether it is inherent 
to galaxies in all environments. Since we see structural bi- 
modality in the most isolated galaxies of the UMa cluster as 
well as in the denser environment of Virgo, it appears that 
this effect is independent of environment. New data from 
field NIR surveys such as PanSTARRS and UKIDSS (e.g.. 
Smith et al. 2008) will be needed to revisit this question for 
field galaxies (provided careful assessment of completeness 
for field samples). 

We now entertain possible physical interpretations for 
structural bimodality. 

5.3 Differences Between Cluster and Field 
Galaxies 

Dressier (1980) first demonstrated the existence of strong 
correlations between various galaxy properties and environ- 
ment. Red sequence galaxies live preferentially in dense en- 
vironments, whilst blue, star-forming galaxies are predomi- 
nantly found in isolation. This dependence on environment 
was reinforced as a major cause for the observed SDSS color 
and star formation bimodality. Not surprisingly, the local 
density of galaxies causes a strong morphological bias. Merg- 
ing and tidal interactions, which are more common in dense 
environments (Moss & Whittle 1993; Gnedin 2003), can de- 
plete the gas content of spiral and irregular galaxies, turning 
them into ellipticals and possibly SOs (Toomre & Toomre 
1972; Farouki & Shapiro 1981; Barnes 1999). An intraclus- 
ter gas can also strip gas from late- type galaxies, reducing 
star formation drastically (Gunn & Gott 1972). Processes 
such as galaxy-galaxy harassment, tidal stripping, and star- 
vation, which are more efficient in dense environments, will 



all reduce the amount of active star formation. It is then not 
surprising that there should be a correlation between color, 
star formation and environment in this context. 

There are also several properties for the dark halo that 
appear to be dependent on environment. The mass, shape, 
spin and formation time of dark matter halos formed in 
denser environments that become clusters or filaments ap- 
pear to be significantly different than those formed in iso- 
lation (Hahn et al. 2007). Simulations suggest that galaxy 
halos with masses M < M,, where M* is the galaxy halo 
mass at which the mass function turns over, have parameters 
that correlate strongly with environment, while those with 
M > M, do not. Hahn et al. (2007) find that low mass halos 
are older and have higher angular momenta when formed 
in proto-cluster environments, while low mass halos in the 
field form later and have lower angular momenta. The de- 
pendence of halo properties on environment could explain 
part of the observed differences in cluster and field galaxies. 

However, using a sample of 329 nearby cluster and field 
star forming galaxies, Vogt et al. (2004) failed to find any 
correlation between stellar mass-to-light ratios or circular 
velocities with environment. More observational and theo- 
retical studies of environmental effects on galaxies will be 
needed. 



5.4 Theoretical Interpretation 

We have found so far that the gas-poor and the gas-rich UMa 
and Virgo cluster galaxies each exhibit a surface brightness 
bimodality. It is unclear whether this result is unique to 
clusters, since there is no complimentary NIR survey of field 
galaxies that reaches the same depth in surface brightness. 
However, UMa very much resembles the field due to its low 
density (see TV97), whilst Virgo is a much richer cluster, 
suggesting that environment may not play a strong role. 

This surprising result requires an explanation which 
current models of galaxy formation fail to provide. However, 
our results may suggest routes for further investigation and 
allow the rejection of certain hypotheses. We now address 
various explanations for the structural bimodality in light of 
the available data and current theoretical ideas. 

(i) Stellar Populations: Previous studies based on SDSS 
imaging for very large number of galaxies have shown bi- 
modal distributions of color and star formation rate (e.g., 
Strateva et al. 2001), which is largely the result of two galaxy 
populations: blue, star-forming spiral galaxies and red, qui- 
escent spheroidal galaxies. We showed in §4.5 that color is 
uncorrelated with the observed surface brightness bimodal- 
ities of cluster galaxies; stellar populations are thus unlikely 
to be a factor. 

(ii) LSB Formation Enhanced in Clusters: If the struc- 
tural bimodality is only present in clusters, it could re- 
sult from more efficient tidal torquing in dense environ- 
ments (Barnes & Efstathiou 1987; Bett et al. 2007; Hahn 
et al. 2007). If dark matter halos acquire more angular mo- 
mentum in a proto-cluster environment, then LSB galaxies 
(which halt their collapse early due to high angular momen- 
tum) would preferentially be formed in these dense environ- 
ments. In this scenario, higher surface brightness galaxies 
could form from the merging of lower SB galaxies if they 
contain gas that can dissipate energy before forming stars. 
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One would thus expect a "smearing" of the LSB peak to- 
wards HSB, but a discrete gap is still unaccounted for. It 
has been well documented that there is a maximum surface 
brightness that galaxy disks will not exceed (Freeman 1970) . 
If this maximal disk surface density is obeyed, then galax- 
ies which evolve at various rates from low to high surface 
brightnesses via merging may ultimately halt their evolu- 
tion at a disk surface density equal to the Freeman value, 
thus enhancing the HSB peak. This scenario does not pre- 
clude the formation of ISB galaxies via the merging of two 
LSB galaxies or other such events, but simply suggests that 
the ISB case might be a rapid transition stage for galaxies as 
they merge and evolve towards the Freeman limit. Indeed, 
we do not see an absence of ISB galaxies in clusters, but 
rather a relative excess of HSB and LSB galaxies. 

Our group is currently studying the velocity function (and 
thus the distribution of angular momenta) of Virgo cluster 
galaxies which will provide further insight into this problem. 
More studies like that of Hahn et al. (2007) would help elu- 
cidate the role of environment on the distribution of various 
dynamical properties. 

(iii) Specific Stable Radial Configurations: Expanding 
upon the ideas of Mestel (1963), TV97 suggested that galaxy 
disks could settle into one of two stable radial configura- 
tions; disks with a dynamically important baryonic com- 
ponent in the center and those without (also known as 
maximal and minimal disks). They suggested that interme- 
diate surface brightness galaxies are transient phenomena 
in high-density environments where interactions re-arrange 
the baryons. This is supported by the fact that, in UMa, 
most intermediate surface brightness galaxies have signif- 
icant neighbors and thus the bimodality is enhanced for 
isolated systems. However, the lack of ISB galaxies in the 
densest regions of the Virgo cluster suggests the opposite 
result. 

If the structural bimodality is unaffected by environment 
then there may be a more fundamental mechanism prevent- 
ing the stability of ISB galaxies. The ISB regime may cor- 
respond to an unstable dynamical configuration, possibly 
implying a bimodality in the spin parameter. A, for disk 
galaxies. 

(iv) Gas Depletion in Clusters Halts Evolution: The 
amount of gas retained by cluster and field galaxies differ 
due to various factors such as high frequency of tidal inter- 
actions in clusters (Moss & Whittle 1993; Gnedin 2003) and 
the dense intracluster medium which can cause gas strip- 
ping (Gunn & Gott 1972). If ISB galaxies form primarily 
via passive evolution and HSB galaxies form primarily via 
mergers, then the cluster environment would act to prevent 
LSB-to-ISB evolution and enhance LSB-to-HSB evolution. 

Without being able to pinpoint the role that environ- 
ment plays in the observed structural bimodalities, it is diffi- 
cult to provide an explanation for their existence. It is likely 
that no single hypothesis above fully explains the surface 
brightness bimodality and that the separate bimodalities in 
early-type and late-type galaxies will require their own ex- 
planations (although their alignment and equal spacing, as 
seen in Table (2] suggests something fundamentally similar!). 
It is difficult to speculate further without a proper analysis 
of the distribution of /i^ for field galaxies, a complimentary 
dynamical analysis of the Virgo cluster, as well as detailed 



numerical simulations for the effects of external processes, 
such as tidal stripping and mergers, or internal processes, 
such as outflows, on surface brightness profiles of cluster 
and field galaxies. 



6 CONCLUSION 

We have constructed the deepest complete sample of near-IR 
photometry of Virgo cluster galaxies to date and with it, we 
were able to achieve a detailed analysis of the distribution 
of surface brightness profiles and luminosities for 286 Virgo 
cluster galaxies. We also extracted griz photometry from 
the SDSS for comparative analysis with a larger sample of 
742 VCC galaxies. This Virgo cluster study also follows from 
our re-analysis of similar data for the smaller UMa cluster 
(McDonald et al. 2008). Our main findings are: 

• There exists a structural bimodality in each of the 
early and late type galaxy populations of the UMa and 
Virgo clusters. From NIR images of these galaxies, we ob- 
serve surface brightness peaks at ^^=21.7 mag arcsec"^ 
and /ig jj=19.7 mag arcsec"^ for late-type galaxies, and 
fil u =19. 7 mag arcsec"^ and /Xe^^=17.8 mag arcsec"^ for 
early-type galaxies. These three peaks are independent at a 
confidence level >90%. 

• The observed structural bimodalities in each galaxy 
class (early or late) is independent of the observed color 
bimodality. 

• The NIR luminosity function for Virgo has a prominent 
dip at intermediate luminosities which is also seen, to a lesser 
extent, in the optical SDSS data. 

The fact that the dichotomy between HSB and LSB 
galaxies is independent of color and morphology implies a 
possible dynamical connection. We have detected a weak, 
but non-negligible, correlation between Vmax and He, but 
more data will be needed to strengthen any dependence. The 
range of densities from UMa to Virgo also suggests that the 
mechanism for producing HSB and LSB galaxies is indepen- 
dent of environment. Without a sample of true field galaxies 
it would, however, be unwise to speculate about any envi- 
ronmental effect on /i^. Unfortunately, current studies of the 
distribution of surface brightnesses of field galaxies are ei- 
ther hampered by dust extinction in the optical or have not 
achieved deep and complete enough coverage at NIR wave- 
lengths. Thus, most reasonable scenarios for field/cluster 
galaxy evolution are still plausible. 

The most obvious and relevant test for the hypothesis 
of different A-distributions in cluster and field environments, 
as proposed by Hahn et al. (2007), is to assemble a complete 
sample of resolved rotation curves for cluster galaxies. This 
would enable new constraints on the distribution of the spin 
parameter. A, which may be responsible for setting /ie (e.g., 
Dutton et al. 2007). Hahn et al. (2007) predict a distribu- 
tion that is skewed towards high values of A in clusters and 
filaments. A complete sample of resolved dynamical data 
for cluster and field galaxies should allow a full assessment 
of the observed structural bimodality in terms of unbiased 
parameters such as dynamical mass and the slope of the 
rotation curve. 

It has been over ten years since the discovery of the 
surface brightness dichotomy in disk galaxies by TV97. We 
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have now extended that study to a broader environment to 
estabUsh the presence of three separate peaks in the distri- 
bution of effective surface brightnesses in all galaxies with 
>90% confidence, and two peaks in the distribution of disk 
central surface brightnesses in disk galaxies with 95% con- 
fidence. The next few years should see significant improve- 
ments in the size and scope (e.g., range of surface bright- 
ness) of galaxy NIR surveys and in studies into the effect of 
environment on fie- The theoretical interpretation for these 
results however awaits a better determination of A for galax- 
ies in the field and in clusters, and a better interpretation of 
the combined effects of merging and stripping to the overall 
surface brightness of galaxies. 
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APPENDIX A: A. MODEL EXPECTATIONS 

In order to characterize the intrinsic distribution of struc- 
tural parameters such as fie, Ve and C28, free from evo- 
lutionary effects, we have generated a Monte Carlo suite 
of idealized galaxy models using combinations of a large 
number of model parameters. A single galaxy light profile 
can be determined with 8 basic parameters. For the bulge 
and disk, the parameters fie,d, ie,d, l-i'e,b, ^e.b and Sersic n 
are selected randomly from Gaussian distributions whose 
widths match roughly observed galaxy distributions. An- 
other 2 parameters govern the location and brightness of 
spiral arms. The location of the arms varies from 1 to 2 disk 
scale lengths and the peak arm brightnesses ranges from 0.3 
to 0.5 mag arcsec"^ above the disk brightness. Finally, all 
model profiles were convolved with a Gaussian function to 
simulate the effect of atmospheric blur; the seeing FWHM 
is our 8th model parameter. We have also forced 1 in every 
5 profiles to be a pure- Sersic function to refiect the presence 
of systems with no disk. The Gaussian distribution of input 
parameters for our models is shown in the top 6 panels of 
Fig. Ell 

In order to infuse more realism into this experiment, 
two different types of error were added to the model pro- 
files: photon (galaxy + sky) random noise and systematic 
sky subtraction error. We determined the appropriate pho- 
ton noise in the H-band from the distribution of measured 
surface brightness errors for our H-band images (McDon- 
ald et al. 2009). From this distribution, an error function 
was determined to add cosmic error to our model surface 
brightness profiles: 

A(^) = 0.00075xexp ~ ^^'^ ^ +0.014 mag arcsec"^(Al) 

This result is very similar to that produced by Courteau 
(1996) for r-band surface brightness profiles. A typical noise 
value ±A(/i) now is added randomly to each model profile 
at each surface brightness level. 

In order to simulate sky-subtraction errors, a system- 
atic "measurement error" was first determined. This value 
is ±0.003% of the sky brightness at H-band. This random 
sky error can artificially produce truncated (over-subtracted 
sky) or upbending (under-subtracted sky) surface brightness 
profiles, depending on its sign. 

We have generated model galaxies by combining the 8 
model parameters, each selected randomly from a realistic 
range of values, and accounting for random and systematic 
errors. The recovered parameter distribution for the 100,000 
simulated surface brightness profiles can be seen in the lower 
3 panels of Fig. lAll for the total galaxy light. The fie distri- 
bution is nearly Gaussian with a peak slightly brighter than 
the disk fie peak due to the inclusion of substantial bulges 
in most profiles. An HSB tail can be seen in the fie distri- 
bution due to severely bulge-dominated systems. These sys- 
tems have higher fie values since bulges typically have higher 
surface brightnesses than disks. The distribution of the re 
parameter is a combination of bulge and disk components; 
the peak is controlled by the bulge Ve for bulge-dominated 
systems and the tail is due to disk-dominated systems. Fi- 
nally, there is a strong peak in the C28 distribution near 
C28=2.8, the expectation value for pure-disk systems. Pro- 
files with C28< 2.8 arise from cored (n < 1.0) systems or 



galaxies with spiral arms. Profiles with C28> 2.8 or higher 
are bulge-dominated. C28 values cover the full range from 
1.5 to 6.5, in excellent agreement with our data (see tj4.4|l . 

With our realistic model surface brightness profiles, we 
can trace the effects of sky measurement errors. This is im- 
portant since such errors can drastically affect the shape of 
the density profile, causing an artificial truncation or anti- 
truncation at large radii. The distribution of total fie could 
thus be biased as a result. In Fig. IA2I we show the effect 
of increasing the sky measurement error from ±0.001% to 
±0.05%. The net effect of a positive sky error is to pull LSB 
galaxies below the sky noise yielding a faint tail, but the 
bright SB peak never shifts. 

Fig. IA3I shows the variations in fie due to over- and 
under-subtracted sky errors. An under-subtraction error 
(upper half) rarely affects the measured fie by more than 0.2 
mag arcsec"^. An over-subtraction error (lower half) has a 
more drastic effect as it can shift a significant part of the 
outer disk below observable levels. Recall that, while these 
tests about sky errors can be illuminating, the measured 
H-band sky error is typically ±0.005% and the effect of a 
systematic sky error is thus minimal. 

Our simulations confirm that if the input structural pa- 
rameters are distributed normally, we should expect a nearly 
normal distribution of effective surface brightnesses. Under 
no conditions can we generate a distribution of fie that devi- 
ates significantly from the normal function. This result can 
be used as an unbiased comparison against galaxy structural 
parameters extracted from real data. 
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Figure Al. Distribution of input parameters (top 6 panels) and measured parameters (bottom 3 panels) for 100,000 model surface 
briglitness profiles. 
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Figure A2. Histogram of effective surface brightnesses of model profiles. The fraction in the upper right is the simulated sky error. 
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Figure A3. Scatter on fie introduced by various sky estimate errors. 
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APPENDIX B: B. GALAXY DYNAMICS 

The study of the distribution of galaxy structural param- 
eters from light profiles can benefit greatly from the avail- 
ability of dynamical information. For instance, TV97 found 
a clear relationship between a galaxy's disk central surface 
brightness and the inner slope of its rotation curve. In order 
to test this claim, and other correlations between dynam- 
ical and photometric parameters, resolved rotation curves 
would be ideal. Failing any information about the shape of 
the rotation curves, we resort to ID measures of maximum 
rotational velocity. Provided sufficient accuracy, they can, 
at least, enable global comparisons with light profile param- 
eters for the spiral galaxies in our sample. 

We retrieved HI data for 58 VCC galaxies 
from the Extragalactic Distance Database (EDD; 
|http :/ /edd. ifa.hawaii.edu/ ) which overlap with our H- 
band sample. These data have been corrected for aperture 
and turbulent effects following TuUy & Fouque (1985) who 
defined the linewidth parameter Wr. We have also verified 
that the linear (almost 1-to-l) relation between radio and 
optically measured line widths observed by Courteau (1997) 
and others is also observed in this sample. All linewidths 
were deprojected as Vtot ~ 0.5W]j = 0.5Wr sin i and sini is 
given by: 

. 2 . 1 - g^ m1^ 
sm I = (Bl) 

1-90 

Here q is the galaxy's axial ratio and 50 = 0.2 is the 
estimated instrinsic thickness of the disk. 

In order to duplicate Fig. 14 from TV97, we are inter- 
ested in computing the ratio Vt/Vtot- Here K is the maxi- 
mum rotation velocity resulting from a purely baryonic sys- 
tem, while Vtot is the observed rotation speed (which in- 
cludes dark matter) at a fiducial radius. TV97 postulate 
that LSB galaxies tend to be dominated by dark matter all 
the way into the innermost regions (T4 ^ Vtot) while HSB 
galaxies are baryon-dominated within roughly two disk scale 
lengths (K ~ Vtot; e.g., Courteau & Rix 1997; D07). 

To compute V*, we use the following equation: 

K{R)^^^^^xAL/L„ (B2) 

where L{< R) refers to the luminosity within the radius, 
R, and M,/Lh is the stellar mass-to-light ratio computed 
using SDSS colors and Bell et al. (2003). The peak value of 
V*(-R) is the value we call V,. 

Fig. IBll shows the bivariate distribution of V,/Vtot as 
a function of disk central surface brightness, I-Io,h for the 
spiral and irregular galaxies with available linewidths from 
the EDD. 

We do not see the strong dynamical bimodality that is 
apparent in TV97, however there does appear to be a cor- 
relation between V*/Vtot and /^o.Hj ^ ^^11 as hints of a 
bimodality in Hq^h- LSB galaxies tend to be dark matter- 
dominated {Vt/Vtot < 0.5), while HSB galaxies would 
be baryon-dominated {Y,/Vtot > 0.5). The galaxies with 
'Vi./Ytot > 1 tend to have bright bulges or compact nucleii 
which cause a spike in the light-weighted rotation curve at 
short radii. However, with linewidths for only 20% of our 
sample, it is impossible to make these claims with statistical 
significance. To do this analysis correctly, we need to either 
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Figure Bl. Bivariate distribution of V*/Vtot and /i* ^ for 58 
disk galaxies in the Virgo cluster for two different mass-to-light 
ratios. The dashed histogram corresponds to the Bell et al. (2003) 
mass-to-light ratios, while the solid histogram corresponds to a 
constant mass-to-light ratio of 1.0. 



obtain reliable linewidths or deep, resolved, optical rotation 
curves for the full sample. Our team is currently pursuing a 
complete characterization of the velocity function of Virgo 
cluster galaxies in order to address these issues. 



© 2009 RAS, MNRAS 000.[llf27l 



